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Abstract 

We present in this paper the analysis of the measurement of the unknown PMNS 
parameters ^13 and 6 at future LBL facilities performing complete three parameters 
fits, each time fully including in the fit one of the atmospheric and solar oscillation 
parameters within its present (future) error. We show that, due to the presence of 
degeneracies, present uncertainties on 623 and Am23 worsen significantly the preci- 
sion on (^13,(5) at future LBL experiments. Only if a precision on the atmospheric 
parameters at least similar to what expected at T2K-I is reached, then the sensi- 
tivities to ^13 and 6 that have been presented in the literature for many facilities 
(where 623 and Am23 are generally considered as fixed external inputs) can indeed 
be almost recovered. On the other hand, the impact on this measurement of the 
uncertainties on the solar parameters, 612 and Amf2 is already negligible. Our anal- 
ysis has been performed using three reference setups: the SPL Super-Beam and the 
standard low-7 /3-Beam, both aiming toward a Mton Water Cerenkov detector lo- 
cated at L = 130 km; the 50 GeV Neutrino Factory with a 40 kton Magnetized Iron 
Detector to look for the "golden channel" — > with baseline L = 3000 km and 
a 4 kton Emulsion Cloud Chamber to look for the "silver channel" z^e with 
baseline L = 732 km. 



1 Introduction 



The results of atmospheric, solar, accelerator and reactor [1] neutrino experiments show 
that flavour mixing occurs not only in the hadronic sector, as it has been known for 
long, but in the leptonic sector as well. The experimental results point to two very 
distinct mass differences^ , Am^^/ ~ 8.2 x 10~^ eV^ and |Am^j^| ?a 2.5 x 10"^ eV^. 
Only two out of the four parameters of the three-family leptonic mixing matrix U pmns 
[4] are known: 9i2 ~ 32° and 623 ~ 45°. The other two parameters, 9is and 6, are still 
unknown: for the mixing angle ^13 direct searches at reactors [5] and three-family global 
analysis of the experimental data [6,7] give the upper bound ^^13 < 11.5°, whereas for 
the leptonic CP-violating phase 5 we have no informations whatsoever. Two additional 
discrete unknowns are the sign of the atmospheric mass difference and the ^23-octant 

(if ^23 7^45°). 

The full understanding of the leptonic mixing matrix constitutes, together with the 
discrimination of the Dirac/Majorana character and the measure of its absolute mass 
scale, the main neutrino-physics goal for the next decade. However, strong correlations 
between ^13 and 6 [8] and the presence of parametric degeneracies in the (^13, S) param- 
eter space, [9]-[12], make the simultaneous measurement of the two variables extremely 
difficult. Several setups have been proposed to face these problems and perform this 
task, the first option being Super-Beam's (of which T2K [13] is the first approved 
one). New machines have been also proposed, such as the /3-Beam [14] or the Neutrino 
Factory [15]. 

In the literature, however, the simultaneous measurement of ^13 and 6 has been nor- 
mally studied considering the solar and atmospheric mixing parameters as external 
quantities fixed to their best fit values (see for example Ref. [16] and refs. therein; 
see also [17] for some recent papers). This is clearly an approximation that has been 
adopted to get a first insight on the problems related to the (^13, 5) measurement. 
However, the experimental uncertainties on these parameters can in principle affect 
the measurement of the unknowns, and it seems important to perform an analysis that 
goes beyond the two-parameters fits presented in the literature. 

In this paper we therefore study, in a systematic way, the impact that "solar" (i.e. 6*12 
and Amf2) and "atmospheric" (i.e. 6'23 and Am23) parameters uncertainties have on 
the measurement of ^13 and 5 at three of the many proposed setups. By doing this we 
want to catch the characteristic features of the inclusion of external parameters un- 
certainties in a (^^13, 5) measurement. A complete six-dimensional fit^ requires a really 
hard computing effort. The authors of Refs. [18,19] obviate this problem marginalizing 
over all the external parameters and reducing the fit to a two-dimentional one. Our 



^ A third mass difference, Am^^^^ ~ 1 eV^, suggested by the LSND experiment [2], has 
not being confirmed yet [3] and will not be considered in this paper. 

^ To which one could add in principle other variables such as the matter parameter or 
systematic errors, [18,20]. 
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approach, conversely, consists of a series of three-parameters fits (taking $12, Amf2, ^23 
and Am|3 in turn as the third fitting variable) to be compared with standard two- 
parameters fits in ^^13 and d. In this way. we realized that the atmospheric parameters 
are the external inputs whose uncertainties are more important in the reconstruction of 
(^13, 5), and that must be better measured in future experiments. We have also tried to 
compare our results with other methods that have been proposed to deal with external 
parameter uncertainties in the measurement of ^13 and S such as the inclusion of a 
covariance matrix in two-parameters x^'s [9] or the so-called CP-coverage [18]. 

We consider here, as exemplificative setups, three CERN-based facilities: 

• the 4 MWatt SPL Super-Beam [21] and a 7 - 100 /3-Beam [22] both aimed at the 
Frejus tunnel where a 440 kT fiducial volume UNO- like Water Cerenkov detector 
[23] could be located with a L = 130 km baseline. 

• the CERN-based 50 GeV Neutrino Factory (see Ref. [16] and refs. therein), with 
two detectors of different characteristics to take advantage of both the "golden" [8] 
and "silver" [24] channels z/g v^ii^t- The two detectors considered are a 40 kT 
magnetized iron detector [25] located at L = 3000 km and a 4 kT emulsion cloud 
chamber [26] located at L = 732 km in the Gran Sasso tunnel. 

By comparing the results at these three, very different, facilities, we deduce that the 
impact of the atmospheric parameters uncertainties is a common problem that future 
experiments looking for ^^13 and 5 will have to face. Of course, this analysis can be 
done for any of the different setups proposed in the literature and not analyzed here. 
Our intention is mainly to address, in this paper, the problem of how uncertainties in 
the atmospheric and solar parameters affect the measurement of (^^13, 5) at setups that 
have been thoroughly discussed than to present a comprehensive comparison between 
two- and three-parameters fits at all of the facilities proposed in the literature. 

The paper is organized as follows: in Sect. 2 we shortly introduce the three facilities 
and the neutrino-nucleon cross-section; in Sect. 3 we remind the central values and 
the uncertainties of solar and atmospheric parameters; in Sect. 5 we review the para- 
metric degeneracies in the measurement of 9i^, and 5 in appearance and disappearance 
channels; in Sect. 4 we introduce the statistical approach used in the paper; in Sect. 6 
we present our results for the measurement of 6*13 and 5 taking into account the un- 
certainties on solar and atmospheric parameters; in Sect. 7 we show the CP-violation 
discovery potential of the considered facilities taking into account the uncertainties 
on atmospheric parameters; in Sect. 8 we eventually draw our conclusions. In App. A 
we compare our statistical approach with other methods; in App. B we present three- 
parameters fits for the three considered setups for different choices of the input pair 
(^"13,^). 
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Fig. 1. Left: (3 -Beam fluxes at the Frejus location (L = 130 km) [28]; Right: SPL Super-Beam 
fluxes at the Frejus location (130 km baseline) [35]. 

2 The experimental setup 

In this section we describe, briefly, the three facihties that we will use in the following 
and we remind the neutrino-nucleon cross-section used throughout the paper. 



2.1 The (3 -Beam 



The /3-Beam concept was flrst introduced in Ref. [14]. It involves producing a beam of 
/3-unstable heavy ions, accelerating them to some reference energy, and allowing them 
to decay in the straight section of a storage ring, resulting in a very intense neutrino 
beam. The chosen ions are ^He, to produce a pure Ve beam, and ^^Ne, to produce a Ve 
beam. We follow the setup proposed in Ref. [22]: the 7 ratio for the two ions has been 
fixed to 7(^He)/7(^^Ne) = 3/5, in order to have both ions circulating in the storage 
ring at the same time; the 7 value has been fixed to 7i8Ne = 100 (i.e., 76He = 60) to 
tune the neutrino/antineutrino mean energy at the maximum of the Ue — >■ 1^^ oscillation 
probability for the CERN to Frejus baseline. A fiux of 2.9x10^*^ ^He decays/year and 
1.1 X 10^*^ ^^Ne decays/year is assumed. Fig. l(left) shows the /3-Beam neutrino fiuxes 
computed at L = 130 km, keeping mg 7^ [27] and taking into account the three 
different decay modes of ^^Ne [28]. The mean energy of the z/g, i^e beams for this setup 
is 0.23 GeV and 0.37 GeV, respectively. Clearly, energy resolution is very poor at such 
low energy, given the influence of Fermi motion and other nuclear effects. Therefore, 
in the following all the sensitivities are computed for a counting experiment with no 
energy cuts [29] . The /?-Beam is a clean environment to produce electron-type neutrinos: 
the main sources of systematic error are the overall flux normalization (that can be 
controlled with a near detector), the deflnition of the flducial volume of the detector 
and the neutrino-nucleon cross-sections. Alternative /5-Beam proposals can be found 
in Refs. [27,30,31,20]. 
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2.2 The Super-Beam 



A Super-Beam is a conventional neutrino beam with a proton intensity higher than 
that of existing (or under construction) beams such as K2K [32], NuMI [33] and the 
CNGS [34] . With respect to the /3-Beam and the Neutrino Factory, neutrino beams of 
a new design, it has the advantage of a well known technology. On the other hand, the 
flux composition (with as the main component for a 7r+ focusing, plus a small but 
unavoidable admixture of i/^, u^. and u^) limits its sensitivity to — > oscillations. 

We follow the setup proposed in Rcf. [21] as a reference: a 2.2 GeV proton beam of 
4 MWatt power (the SPL), with neutrino fluxes computed in a full simulation of the 
beamline in Ref. [35], assuming a decay tunnel length of 60 m. The corresponding fluxes 
are shown in Fig. 1 (right). Notice that this beam was designed, originally, as the first 
stage of a would-be Neutrino Factory, and it has not been optimized as a facility to 
look for — > on its own. Such an optimization has been presented in Ref. [36]. 
Also in this case, as it was for the /9-beam, the main source of systematic error are the 
poorly known neutrino-nucleon cross-sections, the definition of the fiducial volume in 
the far detector and the overall normalization of the flux (with the additional problem 
of new background coming from neutrino species not present in the /3-Beam flux) . For 
this setup, also, we consider two tentative values of systematic error: an "optimistic" 
2% and a "pessimistic" 5%. 



2.3 The Neutrino Factory 



The Neutrino Factory that we consider consists of a SPL-like Super-Beam and a 50 
GeV muon storage ring [37], with 2 x lO^'' muons decaying in the straight section of 
the storage ring per year. Five years of data taking for each muon polarity is envisaged. 
Two detectors of different technology are considered: a 40 kT Magnetized Iron Detector 
(MID) at L = 2810 km; and a 4 kT Emulsion Cloud Chamber (ECC) at L = 732 or 
2810 km. This proposal corresponds to the design of a possible CERN-based Neutrino 
Factory Complex, with detectors located at the Gran Sasso Laboratory (the ECC) and 
at a second site to be defined (the MID and possibly the ECC). Each one of these 
detectors is especially optimized to look for a particular signal: the "golden" channel 
for the 40 kT MID, and the "silver" channel for the 4 kT ECC. The 

corresponding neutrino fluxes are shown in Fig. 2(left). 

The detectors background and systematics for this specific facility have been studied 
in details in Ref. [25] (the Magnetized Iron Detector) and in Ref. [26] (the Emulsion 
Cloud Chamber). 
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Fig. 2. Left: 50 GeV Neutrino Factory fluxes at the Gran Sasso location (L = 732 km) [8]; 
Right: uN and vN cross- sections on water [39]. 

2.4 The neutrino cross-section 



An important source of systematic error is our present poor knowledge of the uN and 
cross-sections for energies below 1 GeV [38]: either there are very few data (the 
case of neutrinos) or there are no data at all (the case of antineutrinos). On top of 
that, the few available data have generally not been taken on the target used in the 
experiments (either water, iron or lead), and the extrapolation from different nuclei is 
complicated by nuclear effects that at the considered energies play an important role. 
For definiteness we show in Fig. 2 (right) the cross-sections on water used for the Water 
Cerenkov detector throughout the paper [39]. Notice that we also used cross-sections 
on iron and lead for the MID and the ECC, respectively. 



3 The leptonic mixing parameters 



In Tab. 1 we remind the values of solar and atmospheric sector parameters used in 
the paper, their present uncertainties and the errors expected after a round of new 
experiments. 

In particular, in the second column of Tab. 1 we report the input values for 612, 623, Am^2 
and Am23 used in the paper. They correspond to the present best fit values for solar 
and atmospheric parameters [40] with the only exception of 623, for which we do not 
use the present best fit value, 623 = 45°, but 623 = 40° to make manifest the im- 
pact of possible octant degeneracies on the results [11]. Notice that throughout the 
paper the experimentally measured atmospheric mass difference (whose present best 
fit value will be labelled as Am^^^) will be fitted with the three-family parameter 
l^'^isl ~ I'^a ~ '^il (s^^ [41] fo^^ different convention). For the solar mass differ- 
ence, on the other hand, we can unambiguously identify the three-family parameter 
Amfg = m-l — ml with the experimentally measured quantity, Am^^z- ^^e third col- 
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umn of Tab. 1 we report the present uncertainties on each of the parameters. Finally, in 
the fourth column, wc present the uncertainties on solar and atmospheric parameters 
that are expected to be achieved with ongoing or planned experiments. For an estimate 
of the reduction of solar parameter uncertainties we refer to Ref . [42] . For an estimate 
of the reduction of atmospheric parameter uncertainties we refer to the Letter of Intent 
of the T2K-phase I experiment, [13]. The expected error on | Am^gl for the central value 
|Am23| — ■^"^otm = 2.5 X IQ-^ cV^ is a function of the sign of the atmospheric mass 
difference, something that will not be measured at T2K-1. For this reason we present 
both spreads specifying the chosen hierarchy, using the results of an analysis yet to 
appear, [43]. 

The T2K-1 improved bounds are used to analyse the impact of the expected atmo- 
spheric uncertainties in the measurement of (^13, 5) at the /5-Beam and the Neutrino 
Factory 3 . On the other hand, in Sect. 5.2 it will be shown that the disappearance 
channel at the SPL Super-Beam can improve significantly the present uncertainties 
on the atmospheric parameters. This measure will therefore be combined with the ap- 
pearance channel when analysing the impact of the atmospheric uncertainties in the 
measure of (6'i3, 5) at the Super-Beam. 



Solar Sector 


Central values 


Present [40] 


Expected (KamLand) [42] 


tan^ 9i2 


0.39 


0.30 - 0.54 


0.30 - 0.51 


Am?2 (eV2) 


8.2 X 10-5 


(7.5 - 9.1) X 10-5 


(7.7 - 8.7) X 10-5 


Atmospheric Sector 


Central values 


Present [40] 


Expected (T2K-phaseI) [43] 


tan^ 6*23 


0.7 


0.53 - 2.04 


0.62-0.85 / 1.21 - 1.66 


|Ami3| (eV^) 


2.5 X 10-3 


(1.7-3.5) X 10-3 


Satm = +, (2.42 - 2.61) X 10-3 
W = -, (2.46 -2.64) X 10-3 








Expected (SPL, this paper) 


tan^ 6*23 


0.7 


0.53 - 2.04 


0.53 - 2.04 




2.5 X 10-3 


(1.7-3.5) X 10-3 


Satm = +, (2.30 - 2.75) X 10-3 
Satm = -,(2.40 -2.90) X 10-3 



Table 1 



Central values and allowed ranges for solar and atmospheric parameters. 



An updated detailed computation of the expected errors on atmospheric parameters that 
can be obtained at this facility is lacking, [44,45]. 
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4 Statistical approach 



In this section we describe the statistical approach used in the paper to estimate the 
impact of uncertainties in the atmospheric and solar parameters in the measurement 

of (9^3,5). 

The obvious approach would be to fit the data in {Na + 2) parameters, with Na the 
number of external parameters that are allowed to vary in a given range (e.g., the 
solar and atmospheric parameters plus the matter density). This procedure, however, 
is increasingly time-consuming as the number of parameters to be fitted goes up. For 
this reason, in order to understand how any single parameter affect the measurement, 
we perform three-parameters fits in ^^13, 6 and one of the following parameters in turn: 
9i2, Am^2) ^23 and Aml^, each of them allowed to vary uniformly in the ranges of 
Tab. 1 (no gaussian priors are introduced). The matter density has been considered as 
a fixed quantity throughout the paper [46] . 

To perform the three-parameters fits we have constructed grids for the expected number 
of charged-current events for each facility, each grid in the two unknowns, ^13 and 5, 
plus the two measured parameters in the solar sector {612, Amf2) or in the atmospheric 
sector {$23, Am^^). When studying the impact of solar parameter uncertainties we 
have fixed the atmospheric parameters to 6*23 = 40° and Am^s = 2.5 x 10~^ eV^ and 
computed four different grids, for Satm = il; Soct = ±1. When studying the impact of 
atmospheric parameter uncertainties we have fixed the solar parameters to $12 — 32° 
and Am?2 = 8-2 x 10"^ eV^. In this case, only two grids must be computed, one for 
each value of Satm- the octant-degeneracy need not to be considered as an external 
(discrete) input, since 623 is one of the free parameters in the grid. 

When a three-parameters fit is performed, the other parameters in the grid are fixed 
to the corresponding present best fit value for Ou, Am^g Am23 ^23 = 40° for 
the atmospheric angle (to take into account possible octant and mixed degeneracies, 
that would disappear for maximal mixing). This procedure is used to study the effect 
of one parameter at a time on the (^13, S) measure. The function is: 



^a^(^13, Satm, Soct) - N^/si^lS, ^, X] Satm, Soct) 



(1) 



with X any of the parameters to be fitted in addition to ^13 and 5, ± refers to neutrinos 
or antineutrinos and A*"^^ is the number of charged leptons observed in the detector 
for a t'a(^'a) beam. The error on the sample A^^^ is: 

(SN^)' - + (4 + (4 B%)' , (2) 

where aj^± is the statistical error on N'^q (Gaussian or Poissonian, depending on the 

a/3 ^ 

corresponding statistics), B^q is the sum of beam and detector backgrounds for the 
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considered channel, computed as in Refs. [21,22], and is the total systematic error 
for the considered channel at a given facility. No covariance matrix for the non-fitted 
parameters has been considered. The three-parameters function defines a three- 
dimensional 90% CL contour that is eventually projected onto the (6'i3, 5) plane to 
perform a direct comparison with the standard two-parameters 90% CL contours for 
the considered setups^ [22,27,28]. 

A discussion on the statistical approach chosen and its difference with existing ap- 
proaches [9,18] is mandatory and can be found in App. A. 



5 Pcirameter correlations and degeneracies 



Parameter correlations and degeneracies arise in the determination of ^13 and 5 at future 
neutrino experiments, as it has been studied in many papers [9]-[12]. The problem is due 

to the strong correlation between these two parameters in the appearance transition 
probabilities (z/g u^, and Ug) and in the present (and near future) ignorance 

of two discrete unknowns, the sign of the atmospheric mass difference Am^g and the 
^23-octant, that can be parametrized by the sign variables Satm — sign[Am|3] and 
Soct = sign[tan(2023)] that take the values ±1 for Am23 > 0(< 0) 623 < 45°(> 45°), 
respectively. Solving the systems of equations corresponding to the four distinct choices 

of Satm and Soct- 

-^0/3(^13) ^atrn, ^oct) = ^afii^lS' ^atrn = Satm] -^oct = ^oct) ; (3) 

-^0/3(^135 ^'1 ^atrm ^oct) = ^afsi^l^^ -^atTn = ~Satm: ^oct = Soct) : (4) 

-^0/3(^13) Satm, Soct) = ^adi^l^^ ^5 Satm = Satm, Soct = —Soct) , (5) 

-^0/3(^13; ^'-1 Satm, Soct) = ^afii^^^^ ^'1 ^atm = —Satm, Soct = —Soct) , (6) 

(with defined in the previous section) will result, in general, in the input pair 
(^13, S) plus seven additional solutions (the clones) to form an eightfold degeneracy: 
the intrinsic clone (Eq. 3), the sign clones (Eq. 4), the octant clones (Eq. 5) and the 
mixed clones (Eq. 6). A complete theoretical analysis of the clones location has been 
presented in Ref. [48], where an algorithm to numerically find each clone location in 
the (^13, S) plane as a function of the considered experimental setup and of the input 
parameters has been given. A similar approach can be applied to study the presence 
of degeneracies in the disappearance channels Ug — > and i'^ ^ i'^- 



5.1 Correlation and degeneracies in disappearance 



The Ue disappearance probability does not depend on the CP violating phase 5 and 
on the atmospheric $23 mixing angle. The ^13 measurement is, therefore, not affected 

^ A preliminary result obtained by means of this procedure has been presented in [47] . 



9 



by (^13 — 5) correlations nor by the Soct ambiguity. The Ve niatter oscillation 

probability, expanded at second order in the small parameters ^13 and (Amf2-^/-E') 
reads: 

(7) 

where A23 = ^nil^/2E, A12 = Am?2/2^, A = y^GpNe and B^ = \A^ A23I with ± 
for neutrinos (antineutrinos) , respectively. This formula describes reasonably well the 
behaviour of the transition probability in the energy range covered by the considered 
/3-Beam setup (L ~ 100 km and Eiy ~ 100 MeV) and it illustrates clearly that two 
sources of ambiguities are still present in z/g disappearance, Satm (for large values of ^13, 
i.e. in the "atmospheric" region) and the (^13 — Ou) correlation (for small values of ^13, 
i.e. in the "solar" region). A /3-Beam could in principle improve our present errors on 
the solar parameters through i/g disappearance. We have checked that this is not the 
case for the considered setup: at large ^13 the second term in eq. (7) dominates over 
the last term, that is more sensitive to solar parameters. On the other hand, for small 
9i3 the statistics is too low to improve present uncertainties on 6*12 and Am^2 (remind 
that energy and baseline of the low-7 /3-Beam has not been chosen to fulfill this task, 
and therefore our results are not surprising at all) . Eventually, in Ref . [49] it has been 
shown that if systematic errors cannot be controlled better than at 5%, the /3-Beam 
disappearance channel does not improve the CHOOZ bound on ^13. 

Eq. (7) can be also applied to reactor experiments aiming to a precise measurement of 
^^13 in a "degeneracy- free" environment. For the typical baseline and energy of a reactor 
experiment (e.g., L — 1.05 km and (Ei,) = 4 MeV for the Double-Chooz proposal, 
[50]) we can safely consider antineutrino propagation in vacuum. As a consequence, no 
sensitivity to Satm is expected at these experiments, since — > A23 for A23 ^ A. It 
is very difficult that reactor experiments could test small values of ^13, and thus the 
^13 — 012 correlation (significant only in the "solar" region) can also be neglected. 



5.2 Correlation and degeneracies in disappearance 

A Super-Beam facility can perform an independent measurement of the atmospheric 
parameters via the disappearance channel: these kind of facilities should in prin- 
ciple reduce the error on the atmospheric mass difference to less than 10 % and on 

the atmospheric angle to ~ 10 % [51]. It is thus interesting to study, as for the z/g 
disappearance channel, the presence of parameter correlations and degeneracies. The 
vacuum oscillation probability expanded to the second order in the small parameters 
^13 and (AuL/E) [52] is: 
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) 



(8) 



where J = cos 6*13 sin 26'i2 sin 26*13 sin 26*23- The first term in the first parenthesis is the 
dominant one and is symmetric under ^23 ^ 7r/2 — 6*23. This is indeed the source of our 
present ignorance on Spct- This symmetry is hfted by the other terms, that introduce a 
mild CP-conserving 5-dependence also, albeit through subleading effects very difficult 
to isolate. We present our results for the disappearance channel in the (^^23) ^^23) 
plane: as a consequence, we do not need to specify the 6'23-octant, since the interval 
O23 £ [36°, 55°] is spanned exphcitly. 

Solving the two systems of equations: 



four different solutions are found for ^23 7^ 45°: two solutions from eq. (9), the input 
value ^23 = ^23 and ^23 — ti'/S — ^23, being the second solution not exactly at ^23 = 
tt/2 — 623 due to the small ^23-octant asymmetry; and two more solutions from eq. (10) 
at a different value of |Am23| [49]. In eq. (8) we can see that changing sign to Am^^ the 
second term becomes positive: a change that must be compensated with an increase in 
|Ami3| to give P^^(Am^i^; Sa^) ~ -^/x/x ( I ^^23 1 ) Satm)- The two solutions of eq. (10) 
corresponding to the wrong choice of Satm can be observed in Fig. 3(left), where equal- 
number- of- events (ENE) curves are computed for both Satm = ^atm (solid) and Satm = 
—Satm (dotted) at the considered Super-Beam facility. The two intersections are notably 
off the input pair ^23 = 40°, Am23 ~ ■^'^Lm- expected, the two sign clones are 
located at |Am23| — ^"^atm almost symmetric with respect to ^23 = 45°. The 

shift in the vertical axis is a function of ^13 and 6. If ^23 = 45° only two solutions 
(corresponding to different choices of Satm) are expected. 

We must also stress that such an uncertainty can be enhanced once we take into 
account that ^13 and S are completely unknown [51] (although the impact of this 
last parameter in z/^ disappearance is expected to be rather small). To gain some 
feeling on the precision that can be expected in a i^^ disappearance measurement at the 
SPL Super-Beam facility, we performed a full three-parameters analysis in ^23, l^^^ilgl 
and 6*13 for the input parameters ^23 — 40°, Am^g = 2.5 x 10""^ eV^ and ^13 = 7°. 
We have then projected the 90% three-parameters CL contour onto the (^^23i 1^^231) 
plane. Fig. 3 (right). In the absence of a complete simulation of the systematics and the 
background for the disappearance channel at the SPL Super-Beam [21], we have 
adopted as an estimate of the expected background and efficiency those used in [22,27] 



^/t(^23, Am^t^; Satm) = ^i(^23, \Amls\; Satm) > 
^/t(^23, Am^^„; Satm) = N^n{023, lAm^gl; -Satm) , 



(10) 



(9) 
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and [28] for the Ve appearance channel at the /3-Beam facihty. A 2% systematic 

error has been assumed. The sohd hne refers to the projection of the three-dimensional 
90 % CL contour on the (6'23, lAm^sl) plane for Satm = Satm^ the dotted line to the 
projection of the 90 % CL contour for Satm = —Satm- 




Fig. 3. Left: ENE curves in the (623, \Am2^\) plane at the SPL Super-Beam facility, for 
Satm = Satm (soUd) and Satm = -Satm (dotted) and 023 = 40°, Am^j^ = 2.5x10"^ eV^. Right: 
Projection of the three-parameters 90 % CL contour in (1923, |Am23|, ^isj o,t the Super-Beam 
facility for the input point 623 = 40°, Am^j„ = 2.5 x 10~^ eV^ and 613 = 7° onto the 
(623, |Am23|j plaiT-^- Again, the solid (dotted) line stands for Satm = Satm (satm = —Satm)- 



As expected, the three-parameters fit presents a second allowed region in the parameter 
space at lAm^gl > Am^^„ when the wrong Satm is considered. Notice that, performing 
a three-parameters fit in ^23, Amgs and 5, the difference between the two- and three- 
parameters contours is much smaller. In [49] it has been shown that a larger spread in 
^23 is found for ^23 7^ 45°. We can perform fits with ^23 non-maximal and for different 
input values for G [0, 10°]. The result of such an analysis is that the SPL Super-Beam 
will be able to measure 6'23 in the interval [36°, 55°] and Am^g in [2.3,2.9] x 10"^ eV^ 
for the input pair ^23 = 40°, Am23 = Am^^^. Notice that the expected SPL precision 
on Am23 is comparable with what expected at T2K-I [13]. On the other hand, the 
expected SPL precision on ^23 is much worse than the T2K-I one, a consequence of 
the fact that the considered SPL setup is a counting experiment and it has no energy 
resolution. 
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Fig. 4. Comparison of the projection of three-parameters 90% CL contours onto the (9i3, 5) 
plane (solid lines) with the corresponding two-parameters 90 % CL contours (dashed lines) 
after a 10 years run at the (3-Beam. Different choices of Satm, Soct oltc plotted separately. The 
input parameters (represented by a star) are: ^13 = 2°, 7°, 5 = 45°. Left panel: x = 812; right 
panel: x = Amf2- 

6 Impact of parameter uncertainties on ^^13 and 6 

In this section, we discuss tlie impact of the uncertainties in the solar and atmospheric 
parameters to the simultaneous measurement of ^13 and S at the considered /3-Beam, 
Super-Beam and Neutrino Factory. The three experiments will be discussed separately. 

In all fits we have combined informations from all available channels for both polarities, 

x\e^,,6,x)=Y.xK0i3,S,x), (11) 

i 

where Xi is the three-parameters x^ function defined as in eq. (1) for a given channel and 
polarity. All channels have been taken as independent measurement and no covariance 
matrix has been introduced, following the approach described in Sect. 4. In general, 
a "pessimistic" systematic error, = 5%, has been used in appearance channels. On 
the other hand, a 2% systematic error has been used in disappearance channels. 

6. 1 The solar sector 

We study the effect of present uncertainties on the solar sector parameters in the 
measurement of ^13 and 5 performing two distinct three-parameters fit in ^13, 5 and 612 
(for fixed Am^g) ^^^12 (^^^ fixed 612). The fits have been performed using 10 years 
of /3-Beam running with both polarities. 

The projection of the three-parameters 90% CL contours onto the (^13, 5) plane are 
presented in Fig. 4. In the left panel we have fixed Am^2 — 8-2 x 10~^ eV^ and drawn 
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the projection of the three-dimensional contours for (^13, 5, ^12), for 612 G [29°, 36°]. 
In the right panel we have fixed 612 = 32° and drawn the projection of the three- 
dimensional contours for x'(^^i3,5,Am22), for Arr^ G [7.5,9.1] x 10"^ eV^. In both 
cases, the atmospheric parameters have been fixed to Am23 ~ 2. 5 x 10~^ eV^ and 
023 — 40°. The input values for the two unknowns are 613 = 2°, 7° and 5 — 45°. 
For each panel, the results for the four different choices of the two discrete variables, 
Satm and Soct, are presented separately. Finally, the projection of the three-parameters 
90% CL contours (solid lines) are directly compared with the two-parameters 90% CL 
contours (dashed lines) obtained fixing the solar parameters to their present best fit 
values, 012 = 32°, Am?2 = 8.2 x 10"^ eV^ 

As we can see in both panels, most of the plotted three-parameters contours coincide 
for any practical purpose with the corresponding two-parameters ones, with small de- 
viations easily explained by the different CL in two- and three-parameters x^- As a 
result, we claim that the impact of solar parameter uncertainties on the measurement 
of ^13 and 5 is negligible for ^13 > 2°. This is indeed a consequence of the sub leading 
dependence of the oscillation probability on the solar parameters (see, for ex- 

ample, Refs. [8,52]) for large values of ^13. When ^13 is large, we are in what has been 
called the "atmospheric regime" in Ref . [9] ; only for ^13 below the verge of the /3-Beam 
^13-sensitivity, i.e. for ^13 > 2°, we enter in the so-called "solar regime". Clearly, no 
signal is expected at the /?-Beam in this case: we can thus safely claim that the solar 
parameter uncertainties do not affect significantly the measurement of ^13 and 5 at the 
considered facihty. 

Similar conclusions can be drawn for the SPL Super-Beam and for different values of 5 
and will therefore not be repeated here. For the rest of the paper, the solar parameters 
will be considered as fixed external inputs: 9i2 — 32° and Am^a = 8.2 x 10~^ eV^. 

6.2 The atmospheric sector at the (5-Beam 

As for the solar sector, we study the effect of present uncertainties on the atmospheric 
sector parameters in the measurement of ^13 and 5 performing two distinct three- 
parameters fit in ^13, 5 and 623 (for fixed Am23) or Am|3 (for fixed ^23)- 

The comparison between two- and three-parameters fits is presented in Fig. 5, where the 
projection of the three-parameters 90% CL contours onto the (^^13, 5) plane (solid lines) 
and the two-parameters 90% CL contours (dashed lines) have been plotted separately 
for each possible choice of the two discrete variables, Satm and Soct- In the left panel we 
have fixed Am^^^ = 2.5 x 10~^ eV^ and drawn the projection of the three-dimensional 
contours for X^(^i3, 5, ^23), for ^23 G [36°, 55°] (see Tab. 1). In the right panel we 
have fixed ^23 = 40° and drawn the projection of the three-dimensional contours for 
X^(6'i3,5, Aml^), for Am^3 G [1.7,3.5] x 10"^ eV^ (see Tab. 1). The input values for the 
two unknowns are ^13 = 2°, 7° and 5 = 45°. The two-parameters contours have been 
drawn using fixed values for the atmospheric parameters, ^23 = 40°, Am^^^ = 2.5 x 10~^ 
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Fig. 5. Comparison of the projection of three-parameters 90% CL contours onto the (9i3, 5) 
plane (solid lines) with the corresponding two-parameters 90 % CL contours (dashed lines) 
after a 10 years run at the (3-Beam. Different choices of Satm, Soct oltc plotted separately. The 
input parameters (represented by a star) are: ^13 = 2°, 7°, 5 = 45°. Left panel: x = 823; right 
panel: x = Am-23. 

eV^. We have checked that no significant improvement is observed when the systematic 
error in the appearance channel is set to 2%. 

In Fig. 5 it is manifest the impact of atmospheric parameter uncertainties on the 
measurement of 613 and 6, for both 623 and Am"^^ fits. Both unknowns are measured 
with errors much larger than those expected from two-parameters contours. This must 
be compared with the results of the previous section, where it has been shown that 
solar parameter uncertainties have a negligible impact. 

Consider first the left panel of Fig. 5: the results from a three-parameters fit in 
(6'i3, 6, 623)- Notice that, being 6*23 a fitting variable in the whole range 623 G [36°, 55°], 
the effect of the "octant ambiguity" is automatically taken into account by the three- 
parameters function. For this reason the three-parameters contours labelled as 
"true" and "octant" are identical. Notice also that this is not the case for the two- 
parameters contours, where the choice of Soct is refiected in contours located at different 
values of 6*13 with respect to the input ^13. For 613 = 7° we can see that a large error in 
6*13 is induced by the uncertainty in 6*23 (with A9i3 as large as 4°). This is a consequence 
of the fact that the leading term in the z/g — > oscillation probability is proportional 
to the combination sin^(2^i3) sin^6'23: to compensate a change in 6'23, a change in 6'i3 
is needed. For smaller values of ^13 this effect is much smaller. The spread in S is, on 
the other hand, extremely similar in two- and three-parameters contours. This is a 
consequence of the fact that a 5-dependence would be induced in the fit through the 
subleading term in the oscillation probability, that is proportional to sin(2^23) and thus 
less sensitive to changes in 6'23 in the almost symmetric interval considered. Notice that 
three-parameters contours have a box-like shape, with no strong 613 — 6 correlation. 
Finally, the largest values of 6*13 are observed for both choices of Satm at lower values 
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of 023, whereas the smallest values of ^13 are reached for larger values of ^23- 

Consider now the right panel of Fig. 5: the results from a three-parameters fit in 
(^13, 5, Aml^). In this case to different choices of Soct correspond different contours (^23 
is a fixed external input and not a free parameter in the fit). As in the previous case, 
a large error in ^13 is induced by the error on the atmospheric parameter, especially 
for ^13 = 7°. The largest value of 6*13 is associated to the smallest value of Aml^, in 
all plots. A characteristic feature of this fit is the significant 5-dependence that can be 
observed in all plots and that was not present in the fits in ^23- The three-parameters 
90% CL contours have a triangular shape (the error in 6 reduces for large values of 
^13), pointing to a strong ^^13 — 6 correlation. It is important to stress that, when Am23 
is the free parameter, the overall error in S is significantly larger in the three- than 
in the two-parameters contours. For both values of ^13 — 2°, 7° roughly half of the 
5-parameter space is covered. This result, considerably worse than what expected from 
two-parameters fit, can be compared with the CP-coverage expectation (explained in 
App. A) for this particular input pair. Fig. 11. For both methods, a rather large error 
in 6 is indeed expected. 

It is clear from Fig. 5 that both atmospheric parameter uncertainties are extremely im- 
portant in the measurement of Oio, and 6: the three-parameters 90% CL allowed regions 
are considerably worse than those obtained with two-parameters fits. In particular, for 
the shown input pairs, 5 would remain completely unknown in the interval S e [0, tt] 
and ^13 would be known only with a large error. 

As an example of how the situation can be improved when using reduced uncertainties 
on the atmospheric parameters, in Fig. 13 we present the projection of the three- 
dimensional 90% CL contours onto the {613,6) plane using the expected uncertainties 
on the atmospheric parameters after T2K-I (last column of Tab. 1): ^^23 £ [38°, 43°] — 
[48°, 52°] and Amig G [2.42, 2.61] x lO'^ eV^ for Satm = + and Aml^ G [2.46, 2.64] x 10"^ 
eV^ for Satm — — , [43]. The octant-ambiguity (that will not be solved at T2K-I) recover 
its discrete nature: separate regions of the parameter space will be spanned by different 
choices of Soct- In this case, all choices of the two discrete variables Satm and Sgct are 
presented together and no comparison with two-parameters contours is shown. In top 
panels x = 6*23; in bottom panels x = Am^g. The results of the three-parameters fit 
with expected uncertainties (right panels), are directly compared with the results pre- 
sented in Fig. 5 computed with the present uncertainties (left panels). The reduction 
of the uncertainties on the atmospheric parameters has indeed an important effect on 
the measurement of ^13 and 6. As it can be seen in the right panels of Fig. 13 a sig- 
nificant reduction of the 6'13-spread is achieved, with plots resembling those obtained 
with standard two-parameters contours and fixed external atmospheric parameters (see 
Refs. [49,28]). The (5-spread is also reduced considerably with respect to the results ob- 
tained with present uncertainties. These comments apply to both ^13 = 2°, 7°. Notice 
that, as expected being ^^23 restricted to one octant only, the octant- and mixed- ambi- 
guities show themselves as separate contours in the {613, 6) plane, as for two-parameters 
fits. 
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A final comment on the impact of the uncertainties on the atmospheric parameters on 
the measurement of ^13 and 5 at the low-gamma /9-Beam is in order. We have shown 
that, with present uncertainties, the measurement of the two unknowns in the PMNS 
mixing matrix is severely spoiled. Errors as large as A6'i3 — 4° are found, and half 
of the parameter space in 5 is spanned for different values of 5. This corresponds to 
a CP-coverage ^ ~ 0.5, a value that spoils completely the possibihty to distinguish a 
CP-violating signal from a CP-conserving one at the considered facility (sec App. A). A 
significant reduction in the uncertainties on the atmospheric parameters is mandatory 
if we plan to use such a facility to look for 5. If 6*23 and ^m^^ can be measured at the 
T2K-I experiment with the expected precision and for any value of ^23, only then the 
results of present two-parameters studies [22,28,49,27] for facilities of this kind can be 
considered rehable. 

In App. B we present the results for different choices of 5, to illustrate the generahty 
of the results above. 

6.3 The atmospheric sector at the Super-Beam 

We now repeat the analysis of the impact of present and expected uncertainties on the 
atmospheric sector parameters in the measurement of ^^13 and 5 at a different facility: 
the SPL Super-Beam. Again, two distinct three-parameters fits in 6*13, 5 and 6*23 (for 
fixed Am23) ^"^23 (^^^ fixed ^23) have been performed, with the Super-Beam running 
2 years with tt"*" and 8 years with tt" to accumulate comparable statistics for neutrinos 
and antineutrinos. A significant difference between this facility and the low-gamma f3- 
Beam considered previously is in that the z/^ disappearance channel at the Super-Beam 
reduces the uncertainties on the atmospheric parameters, as it can be seen in Fig. 3 
(whereas the disappearance channel is useless to this purpose, see Ref. [49]). In this 
case we therefore do not present results using "present" and "expected" uncertainties, 
but we just combine the results from the appearance and disappearance channel. We 
have checked that using a "pessimistic" 5% systematic error in the disappearance 
channel does not change significantly our results. 

The comparison between two- and three-parameters fits is presented in Fig. 6. In the 
left panel we have fixed Am^^^ = 2.5 x 10"'^ eV^ and drawn the projection of the 
three-dimensional contours for X^(^i3) ^23), for ^23 £ [36°, 55°] (see Tab. 1). In the 
right panel we have fixed ^23 = 40° and drawn the projection of the three-dimensional 
contours for x^i^is, Am^g), for Am^g G [2.3, 2.9] x lO'^ eV^ (see Tab. 1 and Sect. 5.2). 
The input values for the two unknowns are ^13 = 2°, 7° and S = 45°. 

The main difference between two- and three-parameters contours resides in that in the 
latter we observe some clones absent in the two-parameters plots. This is a consequence 
of the not satisfactory expected improvement on the error in ^23 and Am^g for ^23 = 40°. 

Also in this case the measurement of 6'i3 and 6 is severely affected by the uncertain- 
ties on the atmospheric parameters. Somewhat smaller errors are found in ^13 and 5 
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Fig. 6. Comparison of the projection of three-parameters 90% CL contours onto the (9i3, 5) 
plane (solid lines) with the corresponding two-parameters 90 % CL contours (dashed lines) 
after a 2-1-8 years run at the Super-Beam. Different choices of Satm, Soct cLfc plotted separately. 
The input parameters (represented by a star) are: 613 = 2°, 7°, 5 = 45°. Left panel: x = 623; 
right panel: x = Am23. 

with respect to the /?-Beam case, but still almost half of the parameter space in 5 is 
spanned for different values of 5. A crucial point is that it does not seem that the 
disappearance channel is capable of a significant reduction in the error on the atmo- 
spheric mixing angle 623. The T2K-I experiment will therefore be crucial, if indeed the 
expected precision in the atmospheric angle can be met for any value of 623. 

In App. B we present the results for different choices of 5, to illustrate the generality 
of the results above. 



6.4 The atmospheric sector at the Neutrino Factory 



We repeat the analysis of the impact of atmospheric parameters uncertainties in the 
measurement of ^13 and 5 at a third facility: the CERN-based SPL-fuelled 50 GeV 
Neutrino Factory. We want to show in this way how the results of Sects. 6.2 and 6.3 
are quite general and must be taken into account at any facility that is considered 
when looking for 6*13 and 5. 

As before, two distinct three-parameters fit in ^13, 5 and 623 (for fixed ^m\^ or /S.m\^ 
(for fixed ^23) have been performed, with the Neutrino Factory running 5 years with 
/i+ and 5 years with 

In the absence of an updated analysis of the expected reduction of atmospheric pa- 
rameters uncertainties at this facility through —>■ z/^, z/g —>■ z/g and z/^j —>■ Vj. (see 
Ref. [44,45] for old analyses), we only present results combining the two appearance 
channels z/g — > z/^ (i.e. the "golden" channel) and z/g — > Vt (i-e. the "silver" channel) 
for both polarities. We use the expected uncertainties after T2K-I, in order to get a 
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Fig. 7. Comparison of the projection of three-parameters 90% CL contours onto the (0i3, 5) 
plane (solid lines) with the corresponding two-parameters 90 % CL contours (dashed lines) 
after a 5+5 years run at the Neutrino Factory. Different choices of Soct o-^^ plotted separately. 
The input parameters (represented by a star) are: 613 = 2°,7° , 6 = 42° . Left panel: x = 623; 
right panel: x = Am23. 

preliminar understanding of the impact of atmospheric parameter uncertainties at this 
facihty. 

The comparison between two- and three-parameters fits is presented in Fig. 7. In the 
left panel we have fixed Am^^^ = 2.5 x 10^^ eV^ and drawn the projection of the three- 
dimensional contours for x^i^is, 5, O23), for 623 G [38°, 43°]-[48°, 52°] (see Tab. 1). In the 
right panel we have fixed 623 = 40° and drawn the projection of the three-dimensional 
contours for x^(^i3, Am^g), for Am^g e [2.4^ 2.7] x IQ-^ eV^ (see Tab. 1 and Sect. 5.2). 
The input values for the two unknowns are ^13 = 2°, 7° and S = 42°. 

First of all notice that at the Neutrino Factory the sign and mixed degeneracies are 
solved, being the magnetized iron detector with a L = 3000 km baseline extremely 
sensitive to matter effects and thus capable to measure Satm- For this reason we only 
present two panels, corresponding to two possible choices of the 6'23-octant, Soct = ^Soct- 
As for the SPL Super-Beam, for small 613 the two- and three-parameters contours prac- 
tically coincide. On the other hand, for ^13 large we must make a distinction between 
Soct = Soct and Soct = —Soct'- whereas the impact of the atmospheric uncertainties for 
Soct = Soct is marginal (something already observed in [9], where the covariance matrix 
approach was adopted and only the right choice of the 6'23-octant was considered), we 
notice how extra octant clones are present in the three-parameters contours that are 
absent in the two-parameters ones when the wrong choice of Soct is taken. This happens 
because in the three-dimensional parameter space O23 cooperates with 6*13 to identify a 
low region with but with ^13 < ^13. 



As for the other facilities, we have seen that the impact of the uncertainties on the 
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atmospheric parameters on the measurement of ^13 and 5 at the Neutrino Factory 
is relevant (albeit perhaps not as important as for the /3-Beam and the Super-Beam 
previously discussed). Again, we stress that the loss in precision is more important for 
large ^13 than for small ^13, a region of the PMNS parameter space that will be selected 
or excluded by the approaching T2K-I experiment. This is indeed a crucial problem 
for precision measurements of the PMNS matrix elements. 

In App. B we present the results for different choices of S, to illustrate the generahty 
of the results above. 



7 CP-violation discovery potential 

Eventually, in Figs. 8-10 we compare the sensitivity to {9i3,6) obtained with a two- 
parameters fit in {9i3,S) or a three-parameters fit in (6*13, 5, 6*23) or {613,6, Am^^) at 
the three considered facihties. The 3a contours have been computed as in Ref . [49] : at 
a fixed ^13, we look for the smallest (largest) value of \S\ for which the two- (three-) 
parameters 3o" contours of any of the degenerate solutions (true, sign, octant and mixed) 
do not touch 5 = 0° nor d = 180°. Notice that, although the input 613 value is fixed, 
the clones can touch 6 = 0°, 180° at ^^13 7^ ^13, also'^ . The outcome of this procedure 
is finally plotted, representing the region in the (^13, 5) parameter space for which a 
CP-violating signal is observed at 3a. Within this approach we can thus take fully 
into account the impact of the parameter degeneracies in the CP-violation discovery 
potential of the three facilities. As for the previous section, we have applied a 2% 
systematic error on disappearance channels and a 5% systematic error on appearance 
channels. As in the previous section we used the expected errors on the atmospheric 
parameters after T2K-1 for the three-parameters fits at the /3-Beam and the Neutrino 
Factory. The SPL Super-Beam analysis relies on SPL data, only (see Sect. 5.2). 

Notice that results are given for the whole allowed range in 5, 5 e [—180°, 180°]. This is 
particularly appropriate, since only an approximate symmetry is observed for \6\ > it/2 
and \6\ < 7t/2 and no symmetry at all between positive and negative 6 in the case of 
the /3-Beam and of the Neutrino Factory. 

Consider first Fig. 8, that refers to /3-Beam results. Notice that the discovery potential 
is not symmetric for positive and negative values oi6, as it has already been observed in 
Ref. [49] . This asymmetric behaviour of the /3-Beam is indeed a statistical mirage caused 
by the low background in the appearance antineutrino sample and the high background 
in the appearance neutrino one (see Ref. [28]). A proper statistical treatment should 
be performed, following Ref. [53], to get rid of this asymmetry for small sin^6'i3: the 

^ This is not the case of Fig. 11 in Ref. [28], where the excluded region in S at fixed ^13 
in the absence of a CP-violating signal at 90% CL is presented. In practice, in that figure 
we compare A^±(^i3, 6) with N±[6i3, 0°), thus obtaining a one-parameter sensitivity plot in 5 
only. 
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Fig. 8. CP-violation discovery potential after 10 years at the low-'j (5-Beam. 
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Fig. 9. CP-violation discovery potential after 2-1-8 years at the SPL Super-Beam. 

treatment, however, is extremely time-consuming and we do not consider meaningful 
applying it here. A further asymmetry can be observed in the different behaviour of the 
three-parameters 3cr contours projection onto the (6'i3, 5) plane for positive and negative 
values of 5: whereas for 5 > we observe that the smallest value of sin^ 6'i3 for which 
a CP- violating phase can be distinguished from a null result goes from [sin^ ^13] mm = 



3 X 10"^ ^ 5 X 10-^ for 5 < we get [sin^ 6 



13 I min 



2 X 10- 



4 X 10^^ for both 



the ^23 and |Am23| fits. A small loss in the discovery potential of this facility with 
respect to the two-parameters fit is observed in both three-parameters fits for negative 
5. In particular, the region in which a CP-violating signal can be distinguished from a 
CP-conserving one goes from 5 e [-80°, -130°] ^ 5 G [-90°, -120°]. 

Consider now Fig. 9, that refers to Super-Beam results. The strong asymmetry for 
positive and negative 5 is not observed in this case, both for two- and three-parameters 
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Fig. 10. CP-violation discovery potential after 5+5 years at the 50 GeV Neutrino Factory. 



fits. The impact of tfie tliird fitting variable, being it 9 



23 



or Am^g, 



is a ratlier small 



loss in the minimum value of sin 6'i3 for which a CP-violating phase is distinguished 

9 X 10^^ 1.2 X 10^3 for ^ ^ Q^,^ and [sin^ ^i3]™„ = 
23. The loss in the (5-interval that is distinguishable 

23 and negligible when 



from a null result: [sin 6*13] mm 
9 X 10"^ 1.4 X 10-3 for x = Amj.^ 

from a null result is rather small for three-parameters fits in Am? 
fitting in 6*23. 



For the Neutrino Factory, Fig. 10, we observe a mixed situation: a strong asymmetry 
between positive and negative S regions (as for the /3-Beam), but a very small difference 
between two- and three-parameters fits (as for the Super-Beam). The asymmetry, how- 
ever, it is not a consequence of asymmetric signal-to-noise ratios^ as for the (3-Bea.m 
but, rather, of a "parametric conspiracy" that for the chosen values of energy and base- 
line results in clones that for many negative values of 6 move toward 5 = 0° or 5 = 180° 
[48], thus preventing a clean identification of a CP- violating signal (see Figs. 22-26). 
The impact of the third fitting variable, being it 623 or Am23, ^ rather small loss in 
the minimum value of sin^ 613 for which a CP- violating phase is distinguished from a 
null result for negative 6: [sin^ ^i3]mm = 1-5 x 10^^ — 2 x 10"^. On the other hand, for 



positive values of 6, [sin 6 



13 Jmm 



2.5 X 10 both for two- and three-parameters fits. 



8 Conclusions 



The simultaneous measurement of 613 and 6 has been often performed in the literature 
considering the solar and atmospheric PMNS parameters as external quantities fixed to 
their best fit values. This is an approximation that has been adopted to get a first insight 

^ It must also be reminded that for 5 years of data taking in each polarity, a smaller statistics 
is accummulated in the wrong-sign muon sample for initial iJ,~ than for initial , due to the 
different i^N and v'N cross-sections. This reduces the sensitivity to ^13 for negative 6. 
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on the problems related to the (^13, 5) measurement. The experimental uncertainties 
on these parameters can in principle affect the measurement of the unknowns, and it 
seemed important to us to perform an analysis that could go beyond the two-parameters 
fits presented in the literature. 

In this paper we therefore have tried to study the impact that solar and atmospheric 
sector parameter uncertainties have on the measurement of ^13 and S at three out of 
the many proposed setups, the standard low-7 /3-Beam, the 4 MWatt SPL Super-Beam 
and the 50 GeV SPL-fuelled Neutrino Factory. By doing this we wanted to catch the 
characteristic features of the inclusion of external parameters uncertainties in a (^13, S) 
measurement. 

Our first goal has been to identify which of the external parameters affects the most the 
results of two-parameters fits. To do so we have performed a series of three-parameters 
fits in 6*13,5 and one of the other parameters (6*12, Am^2) ^23 and Am^j^) in turn as 
the third fitting variable and compared our results with standard two-parameters fits. 
It turned out that the impact of solar parameters uncertainties on the measurement 
of {0i3,5) is negligible, in practice, whereas present uncertainties on the atmospheric 
parameters are large enough to modify in a significant way the results of two-parameters 
fits. In particular, we have noticed that the main cause of the worsening from two- 
to three-parameters fits are the wide displacements of the so-called clones, parametric 
degeneracies due to multiple solutions of eqs. (3)- (6), as a consequence of small changes 
in the external parameters. These results are general to all the considered facihties. 

We have then focused our attention on how the reduction of the atmospheric parameters 
uncertainties could ameliorate the previous results. To this respect, the three facilities 

we have considered are on different footing. On one side, the z/g disappearance channel 
at the standard low-7 /5-Beam cannot improve on its own the present measurement of 
the atmospheric parameters. This facility, therefore, must rely on other experiments to 
meet its goal on ^13 and 5. Luckily enough, it turns out that the precision on $23 and 
Amig expected at the approved T2K-phase I experiment, if met, would be enough to 
improve our three-parameters fits and reproduce the results of two-parameters fits in 
the literature (that, however, were not so good). On the other hand, we have shown 
that the disappearance channel at the 4 MWatt SPL Super-Beam does improve 
the present errors on the atmospheric parameters. This facility, therefore, should not 
necessarily rely on external inputs. The combination of appearance and disappearance 
data, indeed, improve significantly our three-parameters fits. However, 623 is not mea- 
sured well enough and extra clones are still present in the (^13, 5) that are absent in 
two-parameters contours. Finally, the Neutrino Factory has certainly the potential to 
improve significantly the precision on the atmospheric parameters through and i/^ 
disappearance and the — > Vj- appearance channel, something that we have not stud- 
ied in this paper. Using the errors on 623 and Aml^ expected at the T2K-I experiment 
we have checked that extra clones are present in three-parameters fits that were absent 
in the two-parameters analysis. This is a clear indication of the fact that the problem 
we are addressing is common to all the facilities, not only to the low-7 /9-Beam or the 
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SPL Super-Beam. It is not sufficient to just wait and see, but it must be taken into 
account when envisaging future facilities to look for ^13 and S. 

To include the impact of external parameters uncertainties, other methods than direct 

multi-parameter fits have been proposed in the literature. For this reason, we have 
presented a direct comparison of our three-parameters fit results with the so-called 
CP-coverage introduced in Ref. [18]. We have shown that in both methods a significant 
worsening of two-parameters fits arise as a consequence of the inclusion of errors on the 
external parameters in the fit. Whereas the CP-coverage method, however, can be quite 
useful to condense informations about the CP-sensitivity of a facility irrespectively of 
the specific input pair (^13, 5) considered, direct three-parameters fits offer a detailed 
information for both 612 and 6 for specific points in the parameter space. We believe 
that the two methods are, in some sense, complementary and should be combined to 
get a thorough view of the performance of a specific facility designed to measure the 
(^13, S) pair. To this scope we have presented in App. B the results of a series of three- 
parameters fits for the standard low-7 /3-Beam, the 4 MWatt SPL Super-Beam and 
the 50 GeV Neutrino Factory for different choices of the input parameters. 

Eventually, we have studied the impact of the atmospheric parameters uncertainties in 
the CP-violation discovery potential of the three considered facihties. Our results show 
that the discovery potential at the standard low-7 /3-Beam and at the SPL Super-Beam 

is somewhat reduced for negative values of S when uncertainties on $23 and Arri^o^ are 
taken into account. On the other hand the Neutrino Factory appears less affected by 
the inclusion of external parameter errors. 

In conclusion, we think that this paper shows that present uncertainties on atmospheric 
parameters are indeed too large so that the widely adopted approximation of fixing 623 
and Am|3 to their present best fit values be reliable. A new phase of experiments 
that could improve these uncertainties are needed. The precision that is expected on 
the atmospheric parameters at the T2K-I experiment is shown to be such that three- 
parameters fits could reproduce the results of two-parameters fits presented in the 
literature. This experiment is therefore a crucial step in the way to the measurement 
of the two PMNS unknowns, if the precision goals can indeed be met. 

The same kind of analysis we presented here must be clearly repeated at all of the 
proposed setups, something that goes beyond the scope of this paper but is extremely 
important to establish on solid grounds the quest for ^13 and leptonic CP violation in 
the near future. 
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Appendix A 



In this Appendix we review some of the statistical methods that have been proposed 
in the hterature to take into account uncertainties on the external parameters in the 
measurement of {613, S). 

1) Comparison between naive two- and three-parameters fit 

The obvious difference is in the CL contours that can be drawn in the two cases: 
for two-parameters fit the 90% CL corresponds to Ax^ = 4.61, whereas for three- 
parameters fit is Ax^ — 6.25. As a consequence, when a single minimum is found the 
projection of a three-parameters fit onto the two-dimensional contour is, in general, a 
bit larger. The second, not obvious, difference resides in that in the three-parameters 
fit the three-dimensional manifold automatically allows for a displacement of the clones 
solutions arranging for a lower at the relative minima. This is indeed the case for 
the clones corresponding to wrong choices of Satm and of the ^23-octant (sec [48], also). 
If the clones location moves in the three-dimensional manifold, the resulting projection 
onto the plane can be much larger than the two-parameters contour. This is indeed the 
main result of this paper and is discussed at length in Sect. 6. 

2) Inclusion of a covariance matrix in the two-parameters fit 

A fixed error range for any non-fitted parameter can be taken into account introducing 
a covariance matrix in the function as follows: 



xUio^s, 5) = E { [^^(^13, s) - N0^„ 5)] c-^ [N^{e,,, 5) - Nj(e,„ S)] }^^^ 



;i2) 



where Cij is the covariance matrix, i, j refer to different channels at the same experiment 
or to different experiments and {a} is a given set of external parameters. If the errors 
on the entries i and j of the covariance matrix are statistically independent, C is 



da da 



(13) 



where a{a) is the la error on the parameter a. This procedure, followed in [9] for the 
Neutrino Factory and in [22] for the facilities considered in this paper, reproduces the 
enlargement of the two-parameters CL contours observed from a multi-parameter fit 
projected onto the (^13, 5) plane. However, within this approach, the clones locations are 
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not free to move in the multi-dimensional manifold to arrange for a lower x^: they are 
indeed stucked to the location in the (^13, S) plane that can be computed once known 
the external (fixed) parameters (see [48], again, and [54], Sect. 3.3). The displacement 
of the relative minima is indeed the characteristic feature of the multi-parameter fit, 
where Na external parameters cooperate with ^13 and 5 to locate lower regions 
than those found in a two-parameters fit with fixed external parameters. The resulting 
regions are therefore large than those obtained including the covariance matrix in a 
two-parameters fit. 

3) CP-coverage and marginalization over A^^ external parameters 

A parameter that can be used to compare in a condensed way the capability of different 
setups to measure the CP- violating phase 6 has been proposed in [18]. The CP-coverage 
is defined as follows: 

C{S) = Coverage in 5 = t^C/S* M^) , ( 14) 

ZTT 

is the fraction of the 5-parameter space (i.e., 27r) that is allowed at a given CL as a 
result of a measure when the input parameter is 6. The sum goes over N^igg possible 
allowed regions induced by parameter degeneracies, each of them spanning an interval 
A/ (5) of the parameter space. We take the union of these intervals to take into account 
possible partial overlaps of the ^i{S). The smaller the CP-coverage, the better the 
measurement of 5 at a given experiment. In particular, to distinguish a maximally CP- 
violating signal (i.e., 6 = ±90°) at a given experiment from 6 = 0° or 6 = 180°, the 
CP-coverage must be less than 0.5. 

The definition of the CP-coverage above is, however, incomplete: we must still define 
how the A7(^) intervals are computed and which is the dependence of ^(5) on other 
parameters. Indeed, if all the parameters of the PMNS matrix but S were measured, 
^{d) would be just an involute way to express the expected error of an experiment for a 
certain value of 6. This has been called the "CP-pattern", see Fig. 4(right) in Rcf. [54]. 
If, on the other hand, ^13 is also an unknown parameter, we can think of defining a 
function ^(^13, S) and to plot it as a function of different values of ^13 (called "CP- 
scaling" in Ref . [54] ) for a fixed value of S. In this second case, for the particular value 
5 = 0°, the "CP-scaling" as a function of ^13 is nothing else that the CP-sensitivity. 
For example. Fig. 5 in Ref. [22] or Fig. 11 of Ref. [28] represent the sensitivity to S 
for varying ^13 defined as a onc-paramctcr fit where all mixing parameters have been 
measured but 6 and ^^13. The same idea is presented in Fig. 6 of Ref. [22] and Fig. 13 
of Ref. [27] , where the plot represents the capabihty to distinguish a non- vanishing S 
from 5 — 0° or S — 180° at a given one-parameter CL ^ . 

^ It should be stressed that it is not completely correct to present this "CP discovery poten- 
tial" with one-parameter CL contours: being ^13 a parameter to be fitted together with S, we 
should consider two-parameters CL contours, instead. In this case the CP discovery potential 
can be smaller than when only 6 is left as a free parameter, as a result of the fact that the 
two-parameters X^(^i3) ^) can be lower than x^(^i3) ^) for specific values of 9is ^ ^13. This is 
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To take into account the fact that, in general, the parameters of the mixing matrix are 
known only with a finite precision and that 613 is completely unknown at present (and 
thus a one-parameter 5-sensitivity plot is generally overestimating the performance of a 
given experiment), the authors of Ref. [18] have proposed to compute the CP-coverage 
^((5) by first minimizing a {Na + 2)-parameter over external parameters. In this 
way, for any input pair ^13 and 5, a two-dimensional surface of the minimum as 
a function of ^13 and 5 is generated. If we then minimize the resulting function in 
^13, also, we can deduce a one-dimensional function of 6 and of the input parameters 
representing the minimum value of the {Na + 2)-dimensional for ^ given value of S. 
From this marginalized we can finally compute the allowed A/(5) intervals imposing 
that Xmini^j ^) be equal to a given one-parameter CL. 

Clearly, this procedure can fail when multiple minima are present at each marginal- 
ization step. When multiple minima are present, the minimization procedure will gen- 
erally look for the absolute minimum. In this way, the information on other relative 
minima in the can be lost. This is particularly problematic since we know that, 
due to parametric degeneracies, multiple minima should be present. Marginalization 
near a second minimum will give a second function Xmini^y from which a new set of 
A/ (5) intervals can be found. The procedure suggested in [54] is indeed to marginalize 
around each of the expected minima in the multi-dimensional x^ ^^^^ to draw several 
distinct one-dimensional functions Xmm(^'^)- Imposing on any of them the constraint 
Xmm(^> ^) — the full set of allowed regions in 5 at a given one-parameter CL for a 
fixed input 5 is deduced and the ^{5) can be finally computed. Since the minimization 
procedure must be repeated several times (once per expected minima), it is useful to 
choose the starting point of the minimization algorithm in a clever way. In [54] it is 
suggested to solve eqs. (3)- (6) as it has been done in Ref. [48] to find the expected 
clone locations and to use the latter as starting points for the minimization. Applying 
the previous algorithm it is possible to deduce a ^(5, ^13) parameter that reduce signifi- 
cantly the overestimation of the experiment performances in the measurement of 5 that 
is typical of the one-parameter ^-sensitivity plots. We should therefore compare this 
procedure with the projection onto the 5-axis of our three-parameters fits to under- 
stand if a residual overestimation is still present. To this purpose, in Fig. 11 we present 
the CP-coverage ^(^13, 5) and the projection onto the 5-axis of the three-parameters CL 
contours for different values of ^13 and 5. The results in the figure have been obtained 
using the considered /3-Beam facility, for definiteness. 

As it can be seen in the figure, some underestimation of the error on 5 at the considered 
experiment is still present when computing ^{5) and comparing it with the 5-axis 
projection of the three-parameters fits. The main interest of the CP-coverage parameter 
is that the algorithm described above can be iterated for any value of 613 to obtain a 
"CP-scaling" for any given value of 5, thus replacing the (poorly reliable) one-parameter 
5-sensitivity plots. This will give the general picture of the expected error on 5 at a 
given facility, to be complemented in our opinion with multi-parameter fits to particular 

indeed what reported in Fig. 7 of Ref. [49] . 
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Fig. 11. Comparison of the CP-coverage ^ and the 5-axis projection of the three-parameters 
90% CL contour to the same input parameters at the f}-Beam. The marginalized Xmrnl*^'^) 
function is plotted for different choices of Satm and Soct: Satm = Satm,Soct = Soct 

(solid); Satm = -Satm, Soct = Soct (dotted); Satm = Satm, Soct = -Soct (daskcd); 

Satm = —Satm, Soct = —Soct (dot-dashcd) . Horizontal thick black lines are the A/(5) inter- 
vals and the gray-shaded region are the corresponding three-parameters 90% CL contours. Of 
the two possible three-parameters fit (i.e., in 623 or Am2-^), that with the largest error in 5 
is reported. The thin vertical line displays the value of 5 for each plot. The different plots 
refer to: = 2_°,5 = 0° (top left); 613 = 2°,6 = 45° (top right); 613 = T ,5 = 0° (bottom 
left); = 7°,5 = 45° (bottom right). The atmospheric input parameters are 823 = 40°, 
Am2,^ = 2.5 X 10-3 eV. 

values of 6*13 and 6 to get a robust estimate of the facility performance. 
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Appendix B: three-psirameters fits 



In this Appendix we present the projection of the three-dimensional 90% CL contours 
onto the (^^13, d) plane for the three reference setups (the low-7 /9-Beam, the SPL Super- 
Beam and the 50 GeV Neutrino Factory) for different choices of the input pair {6is, 6). 
The external input parameters in all plots are: 9i2 = 32°, Aml2 = 8.2 x 10~^ eV^; 
O23 — 40°, Amig = 2.5 X 10"^ eV^. In this case, all choices of the two discrete variables 
Satm and Soct ^rc presented together and no comparison with two-parameters contours 
is shown. 

For the /3-Beam plots we compare the impact of the present uncertainties on the at- 
mospheric parameters (third column of Tab. 1) with that of the expected uncertain- 
ties after T2K-I (last column of Tab. 1): ^23 e [38°, 43°] - [48°, 52°] and Aml^ G 
[2.42,2.61] X 10-^ eV^ for Satm = + and Am^g G [2.46,2.64] x 10"^ eV^ for Satm = -, 
[43]. Notice that the error on 623 is just the expected error at T2K-I, [13], shifted 
around 623 = 40°. It has been shown in Sect. 5.2 that the u^^ disappearance channel at 
the Super-Beam is rather effective in reducing the uncertainties on the atmospheric pa- 
rameters (whereas the Ug disappearance channel at the low-7 /?-Beam is useless to this 
purpose, see Ref. [49]). In this case we therefore do not present results using "present" 
and "expected" uncertainties, but we just combine the results from the appearance 
and disappearance channel. Finally, for the Neutrino Factory we have considered only 
the two appearance channels z/g v^, 1/^ with the atmospheric parameters with the 
expected uncertainties after T2K-1. 

In general, a "pessimistic" systematic error, e''^ = 5%, has been used in appearance 
channels. A 2% systematic error has been used in disappearance channels. However, 
we checked that using a "pessimistic" 5% systematic error in the disappearance channel 
does not change significantly our results. 

The input parameters are: = 2°, 7°; 5 = 90°, 45°, 0°, -45°, -90°. 
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Fig. 12. Three-parameters 90 % CL contours after a 10 years run at the (3-Beam. Input 
parameters: Oy^ = 2°, 7°; 5 = 90°. Top panels: x = 623; bottom panels: x = Am^s. Left 
panels: present uncertaintes; right panels: after T2K-I. 
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Fig. 13. Three-parameters 90 % CL contours after a 10 years run at the (3-Beam. Input 
parameters: Oy^ = 2°, 7°; 5 = 45°. Top panels: x = 623; bottom panels: x = Am^s. Left 
panels: present uncertaintes; right panels: after T2K-I. 
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Fig. 14. Three-parameters 90 % CL contours after a 10 years run at the (3 -Beam. Input 

As- 



parameters: 9i3 = 2°,7° ; 6 = 0° . Top panels: x = 623; bottom panels: x = Am^s. Left panels: 



present uncertaintes; right panels: after T2K- I. 
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Fig. 15. Three-parameters 90 % CL contours after a 10 years run at the (3-Beam. Input 
parameters: ^13 = 2°, 7°; 6 = —45°. Top panels: x = 623; bottom panels: x = Am^s. Left 
panels: present uncertaintes; right panels: after T2K-I. 



36 





Fig. 16. Three-parameters 90 % CL contours after a 10 years run at the (3-Beam. Input 
parameters: ^13 = 2°, 7°; 6 = —90°. Top panels: x = 623; bottom panels: x = Am^s. Left 
panels: present uncertaintes; right panels: after T2K-I. 
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Fig. 17. Three-parameters 90 % CL contours after a 2+8 years run at the Super-Beam. Input 
parameters: 613 = 2° ,7° ; 6 = 90°. Left: x = 623; right: x = Am^s. 
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Fig. 18. Three-parameters 90 % CL contours after a 2+8 years run at the Super-Beam. Input 
parameters: ^13 = 2°, 7°; 5 = 45°. Left: x = 623; right: x = Am23. 
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Fig. 19. Three-parameters 90 % CL contours after a 2+8 years run at the Super-Beam. Input 
parameters: 613 = 2°, 7°; 5 = 0°. Left: x = 623; right: x = Am^s- 
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Fig. 20. Three-parameters 90 % CL contours after a 2+8 years run at the Super-Beam. Input 
parameters: 613 = 2°, 7°; 5 = —45°. Left: x = 623; right: x = /^m^^. 
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Fig. 21. Three-parameters 90 % CL contours after a 2+8 years run at the Super-Beam. Input 
parameters: 613 = 2°, 7°; 6 = —90°. Left: x = 623; right: x = Am23. 
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Fig. 22. Three-parameters 90 % CL contours after a 5+5 years run at the Neutrino Factory. 
Input parameters: 613 = 2°, 7°; 8 = 90°. Left: x = 623; right: x 
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Fig. 23. Three-parameters 90 % CL contours after a 5+5 years run at the Neutrino Factory. 
Input parameters: 613 = 2° ,7° ; 6 = 45°. Left: x = 623; right: x = Am23- 
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Fig. 24. Three-parameters 90 % CL contours after a 5+5 years run at the Neutrino Factory. 
Input parameters: ^13 = 2°, 7°; ^ = 0°. Left: x = 623; right: x = Am23. 
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Fig. 25. Three-parameters 90 % CL contours after a 5+5 years run at the Neutrino Factory. 

^23- 



Input parameters: ^13 = 2°,7° ; 6 = —45°. Left: x = 623; right: x = Am23. 
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Fig. 26. Three-parameters 90 % CL contours after a 5+5 years run at the Neutrino Factory. 
Input parameters: ^13 = 2°, 7°; 5 = —90°. Left: x = 623; right: x = Am2^. 
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